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I.  INTRODUCTION  AND  PURPOSE 


To  my  mind,  there  are  two  great  unexplained  mysteries  in  our 
understanding  of  the  universe.  One  is  the  nature  of  a  unified 
general  field  theory  to  explain  both  gravitation  and 
electromagnetism.  The  other  is  an  understanding  of  the  nature  of 
turbulence.  After  I  die,  I  expect  God  to  clarify  general  field  theory 
for  me.  I  have  no  such  hope  for  turbulence.  (Brown  1991) 

This  quotation  rather  elegantly  expresses  the  complex  problem  of 
turbulence  in  fluid  flows,  though  its  origin  is  a  matter  of  considerable  speculation. 
Variously  attributed  to  the  Hungarian-American  aerodynamicist  Theodore  von 
Karman,  the  British  mathematician  Horace  Lamb,  and  the  German  theoretical 
physicist  Werner  Heisenberg,  it  remains  fundamentally  appropriate  some  half 
century  after  its  initial  pronouncement  despite  decades  of  advances  in  both 
quantum  physics  and  fluid  dynamics. 

In  the  classical  application  of  Newton’s  second  law  to  wave  dynamics, 
strictly  nonlinear  interactions  between  perturbations  in  the  flow  are  trivially  small 
compared  to  other  forcing  terms  and  are  usually  neglected  (Stull  1988;  Holton 
2004).  It  is  precisely  these  nonlinear  terms,  however,  that  define  the  essential 
mechanics  of  turbulence  in  fluid  flows.  Given  the  mathematical  difficulty  the 
inclusion  of  higher-order  nonlinearities  imparts  to  the  rigorous  theoretical 
treatment  of  turbulence,  numerical  simulation  of  turbulent  flow  is  thus  dependent 
to  some  extent  upon  decades  of  empirical  studies  consisting  of  both  in  situ 
experimentation  and  computer  modeling  endeavors.  Such  studies  have  yielded 
numerous  techniques  for  parameterizing  the  bulk  exchanges  of  mass  (including 
water  and  solid  constituents),  momentum,  and  energy  between  the  Earth’s 
surface  and  the  atmosphere  (Stewart  1979). 

Although  the  ubiquitous  Navier-Stokes  and  mass  continuity  equations 
theoretically  describe  the  physics  of  turbulence,  its  complex  effects  elude 
expression  in  explicit  detail  due  to  the  inherent  nonlinearities  associated  with 
turbulent  flow.  In  the  particular  case  of  the  Earth’s  atmosphere,  turbulence 
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imposes  apparent  randomness  on  the  fields  of  mass,  momentum,  and  energy; 
and  the  resulting  system  evolves  chaotically  in  time  and  space.  Given  the 
formidable  problem  of  not  only  characterizing  the  initial  state  and  boundary 
conditions  for  turbulence,  but  also  developing  methods  to  iterate  its  evolution,  the 
treatment  of  turbulence  in  numerical  weather  prediction  has  historically  taken  a 
stochastic  approach  (Lenschow  1986).  In  this  thesis,  one  characteristic  of 
turbulence  in  the  lower  atmosphere — the  roughness  length — is  calculated  using 
observations  over  heterogeneous  surfaces  and  correlated  to  the  physical 
properties  of  landscapes  near  the  observation  site.  To  facilitate  this  study,  a  brief 
discussion  of  useful  terminology  and  concepts  follows. 

A.  THE  ATMOSPHERIC  BOUNDARY  LAYER 

The  portion  of  the  lower  atmosphere  extending  vertically  from  the  Earth’s 
surface  to  a  variable  height  between  100  m  and  3  km  comprises  the  atmospheric 
boundary  layer  (ABL)  (Stull  1988);  so  called  because  it  loosely  defines  the  fluid 
boundary  between  the  nongaseous  surface  of  the  planet  and  the  remainder  of 
the  atmosphere  (hereinafter  the  free  atmosphere).  More  precisely,  the  generally 
accepted  qualitative  definition  of  the  ABL  is  that  portion  of  the  troposphere 
recently  in  close  proximity  to  the  Earth’s  surface  whose  internal  flow 
characteristics  respond  to  forcing  from  surface  features  on  time  scales  of  up  to 
one  hour.  In  the  foregoing  sentence,  the  word  recently  implies  “within 
approximately  the  previous  day,”  to  acknowledge  the  dominance  of  the  diurnal 
cycle  on  physical  processes  within  the  ABL  (Stewart  1979). 

The  ABL  contains  several  component  layers  distinguishable  by  their 
physical  properties.  Closest  to  the  surface,  usually  extending  no  higher  than  a 
few  centimeters  in  depth  is  the  interfacial,  or  viscous  layer.  It  is  only  in  this  thin 
envelope  of  air  most  directly  in  contact  with  the  surface  that  transport  of  physical 
fields  via  turbulence  is  less  efficient  than  by  molecular  transport  (Stull  1988). 
Above  the  interfacial  layer,  the  lowest  10%  to  15%  (by  height)  of  the  ABL  is 
called  the  surface,  or  constant  flux  layer.  In  this  layer,  fluxes  of  momentum  and 
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energy  are  approximately  constant  in  the  vertical,  and  variations  in  the  wind  field 
are  most  directly  attributable  to  static  stability  and  local  frictional  effects.  The 
influence  of  the  Earth’s  rotation  is  insignificant  and  often  neglected  in  the  surface 
layer  (Kaimal  and  Finnigan  1994). 

The  outer,  or  Ekman  layer  comprises  the  bulk  (about  90%  by  depth)  of  the 
ABL.  Flow  in  the  outer  layer  usually  exhibits  less  variability  in  both  direction  and 
speed  than  in  the  surface  layer,  hence  the  Coriolis  effect  becomes  a  necessary 
consideration  (Stull  1988).  Moreover,  in  the  outer  layer,  fluxes  of  energy  and 
momentum  may  vary  significantly  with  height.  If  sufficient  lift  and  water  vapor  is 
present,  cumuliform  clouds  may  develop  in  this  layer,  which  alter  the 
thermodynamic  response  of  the  entire  ABL  (Wang  2009).  The  entrainment  zone 
divides  the  ABL  from  the  free  atmosphere.  With  few  exceptions,  this  layer 
defines  the  upper  boundary  for  turbulence  generated  in  the  ABL.  The  high  static 
stability  of  this  layer  often  results  in  a  thermal  inversion,  characterized  by  an 
increase  in  temperature  with  height. 

B.  ABL  PARAMETERIZATION  IN  NUMERICAL  MODELS 

1.  Turbulent  Transport 

Recent  advances  in  supercomputing  technology  have  ushered  in  an  era  of 
atmospheric  modeling  on  ever-smaller  domains  and  at  finer  spatial  and  temporal 
resolutions.  It  has  become  routine  practice  for  researchers  studying  microscale 
phenomena  such  as  sea  breezes  and  convective  storms  to  employ  dynamical 
models  with  grid  points  spaced  less  than  1  km  (Kalnay  2003;  Rao  and  Fuelberg 
2000).  Such  scales  are  dimensionally  comparable  to  the  largest  turbulent 
structures  in  the  lower  atmosphere;  therefore,  selecting  optimal  methods  of 
characterizing  the  complex  thermodynamic  adjustments  caused  by  turbulence 
becomes  increasingly  important  at  these  resolutions.  It  is  also  readily  apparent 
that  turbulent  transport  of  fields  of  meteorological  interest  occurs  on  length  scales 
that  are  much  smaller  than  1  km  (Kalnay  2003).  To  account  for  these  subgrid- 
scale  processes,  modelers  usually  parameterize  their  net  effect  at  the  resolution 
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of  the  model.  Researchers  have  developed  numerous  schemes  for 
parameterizing  subgrid-scale  processes  in  the  atmosphere,  and  this  is  an  area  of 
ongoing  research.  Ideally,  as  our  understanding  of  processes  associated  with 
atmospheric  turbulence  continues  to  improve,  so  will  our  skill  in  accurately 
simulating  these  processes  in  numerical  models. 

2.  Roughness  Length 

a.  Definition  and  Measurement  Methods 

Aerodynamic  roughness  length,  frequently  denoted  in  the  literature 
by  z0 ,  is  the  height  above  the  local  solid  or  liquid  surface  at  which  the  velocity  of 

the  mean  wind  is  zero  (Stull  1988).  By  this  definition,  roughness  length,  although 
defined  at  a  specific  point,  substantially  depends  upon  the  morphology  of  the 
surface  in  some  neighborhood  surrounding  that  point.  Surfaces  with  taller, 
irregularly  spaced  individual  roughness  elements  (trees,  buildings,  etc.)  generally 
induce  larger  turbulent  structures  and  thus  have  a  larger  z0,  while  relatively  flat 

surfaces  scarcely  disturb  the  flow  and  correspond  to  smaller  z0.  Through 

decades  of  experimental  work,  over  various  types  of  land  surfaces,  ABL 
researchers  established  tables  of  roughness  lengths  arranged  according  to  land 
use  categories  (Figure  1 ). 
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According  to  de  Rooy  and  Kok  (2004),  the  influence  of  variable 
roughness  reaches  its  maximum  at  the  height  of  local  z0  and  vanishes  at  the  so- 

called  blending  height  (Wieringa  1976,  1986).  The  experimental  work  of 
Wieringa  (1986)  and  Caton  (1977)  suggests  that  over  various  surfaces,  the 
influence  of  roughness  elements  is  appreciable  up  to  about  3  km  downstream 
and  over  a  20°  to  30°  sector.  In  addition,  Wieringa  points  out  that  noticeable 
variation  in  downstream  gustiness  often  occurs  on  a  seasonal  cycle,  presumably 
a  consequence  of  variable  foliation. 
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Roughness  over  the  open  ocean  and  large  lakes  varies  with  wave  height. 
Charnock  (1955)  established  the  empirical  expression  that  relates  surface  wind 
stress  to  wave-induced  roughness.  As  wind  stress  increases,  ocean  waves  grow 
higher,  consequently  roughness  length  increases.  In  this  expression, 


zo  =  a — ,  (1) 

g 

ut  (introduced  in  Chapter  II)  and  z0  have  their  customary  interpretation,  g  is  the 

acceleration  due  to  gravity  at  the  surface  of  the  ocean,  and  the  dimensionless  a  , 
known  as  the  Charnock  parameter,  has  the  experimentally  determined  value  a  » 
0.0144,  though  research  suggests  a  also  has  a  dependence  on  sea  state 
(Komen  et  al.  1998). 

b.  Roughness  Length  in  Numerical  Models 

Successful  prediction  of  the  near-surface  wind  field  in  numerical 
models  is  contingent  upon  accurately  characterizing — either  explicitly,  or  more 
commonly,  by  means  of  one  or  more  parameterization  schemes — the  variation  in 
fluxes  of  momentum  between  the  atmosphere  and  the  Earth’s  surface.  Since 
surface  frictional  stress  is  a  primary  sink  for  atmospheric  momentum  (Newton 
1971),  accurate  parameterizations  of  surface  roughness  improve  model 
verification  throughout  the  entire  atmosphere.  Roughness  length,  like  many 
other  subgrid-scale  processes,  is  typically  averaged  over  some  finite  portion  of 
the  model  domain,  such  as  the  area  contained  in  one  grid  box  (Kalnay  2003). 
However,  observations  of  low-level  dynamic  fields,  particularly  wind,  often  reveal 
strong  dependences  on  local  surface  features  on  scales  far  below  the  resolution 
of  the  model.  Modelers  regard  this  discrepancy  between  observed  and 
forecasted  values,  or  representation  mismatch,  as  part  of  the  inherent  error  of  the 
model.  Jacobs  and  Maat  (2005)  stated  that  for  the  10-m  wind  field,  this  primary 
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source  of  this  error  is  the  difference  between  the  actual  surface  roughness  and 
that  used  in  the  model’s  parameterization  scheme. 

Currently,  numerical  weather  prediction  models  incorporate  surface 
roughness  as  a  fixed  parameter  for  each  grid  point  in  numerical  weather 
prediction  models  (Wang  2009).  Such  a  simplistic  treatment  may  be  appropriate 
when  modeling  the  slowly  varying  distribution  or  characteristics  of  surface 
roughness  elements,  (e.g.,  land  use  changes  or  seasonal  variations  in  leaf 
volume  or  snow  cover).  However,  research  conducted  by  compiling  and 
analyzing  observational  data  presented  herein  suggests  that  such  simplification 
may  be  unsatisfactory  in  all  cases.  Specifically,  if  roughness  length  were 
determined  empirically  at  a  grid  point  near  the  boundary  between  regions  having 
different  roughness  characteristics,  such  as  along  a  coastline  or  at  the  edge  of  an 
urban  area,  it  is  reasonable  to  expect  that  different  values  would  be  derived  if  the 
near-surface  wind  changed  direction.  If  the  upstream  surface  were  smoother,  it 
would  induce  smaller  turbulent  structures  in  the  flow,  which  should  contribute  to 
an  effective  decrease  in  roughness  length  at  the  point  of  measurement.  If  at  the 
same  point,  the  wind  were  to  change  direction  such  that  air  passing  the  point 
originated  over  a  rougher  surface,  it  is  reasonable  to  expect  larger  disturbances 
in  the  wind  field  and  consequently,  a  larger  roughness  length. 

3.  Internal  Boundary  Layers  over  Heterogeneous  Surfaces 

Often,  the  surface  layer  contains  one  or  more  internal  boundary  layers 
(IBLs)  bounded  above  by  a  discontinuity  in  some  variable  or  state  of  the  surface 
layer  (Garratt  1992).  IBLs  result  from  flow  across  surfaces  that  vary  in  any 
quantity  capable  of  affecting  the  physical  properties  of  the  atmosphere  (e.g., 
temperature,  humidity,  or  surface  roughness).  Figure  2  is  a  schematic  diagram 
that  illustrates  the  concept  of  an  IBL  caused  by  a  change  in  surface  roughness. 
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Figure  2.  Diagram  of  an  internal  boundary  layer.  From  Stull  (1988). 

Here,  the  height  of  the  IBL,  8,  is  a  function  of  fetch,  the  distance 
downstream  from  the  point  roughness  changes.  In  cases  of  significant 
heterogeneity,  multiple  IBLs  may  be  present,  each  in  response  to  an  upstream 
surface.  It  is  therefore  important  to  recognize  the  presence  of  IBLs  and  interpret 
observations  made  over  heterogeneous  surfaces  with  care  since  Monin-Obukhov 
similarity  fails  in  the  transition  region  between  an  IBL  and  the  advected  boundary 
layer. 


4.  Surface  Flux  Parameterization  over  Heterogeneous  Surfaces 

Early  research  in  simulating  surface  fluxes  in  general  circulation  models 
(GCMs)  mainly  focused  on  parameterizing  bulk  exchanges  of  mass,  momentum, 
and  energy  between  homogeneous  surfaces  and  the  ABL.  The  resolution  of 
early  GCMs  was  on  the  order  of  tens  to  hundreds  of  kilometers,  and  necessarily 
dealt  with  surface  roughness  heterogeneities  at  the  subgrid  scale. 

Louis  (1979)  devised  a  parameterization  scheme  for  turbulent  fluxes  that 
included  static  stability  effects  and  evaluated  its  performance  in  a  10-day  forecast 
model.  Due  to  its  relative  simplicity  and  brief  calculation  time,  model  developers 
incorporated  this  scheme,  with  occasional  modifications,  into  many  numerical 
models  in  the  years  following.  Wang  et  al.  (2002)  proposed  adaptations  to 
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Louis’s  approach  to  better  characterize  fluxes  over  smooth  surfaces  and  in  cases 
where  the  effective  roughness  lengths  for  momentum  and  heat  fluxes  differ  due 
to  varying  stability. 

Beljaars  and  Holtslag  (1991)  considered  the  effects  of  variations  in 
vegetative  cover  in  aerial  averaging  techniques  for  parameterizing  surface-layer 
fluxes.  Their  work  showed  a  strong  correlation  between  roughness  length  and 
surface  heterogeneities  located  up  to  5  km  upwind  of  their  observation  site  and 
suggested  that  wind  gustiness  could  determine  the  effective  local  roughness 
length  at  such  scales. 

Two  methods  are  commonly  used  to  account  for  subgrid-scale  surface 
heterogeneity  in  GCMs.  The  technique  of  parameter  aggregation  uses  the 
fractional  coverage  of  different  surface  types  within  a  grid  box  to  obtain  grid- 
averaged  parameters  and  then  incorporates  the  subgrid  averages  into  the  bulk 
flux  parameterization  for  the  larger  grid.  The  flux  aggregation  method  uses  bulk 
parameterization  to  obtain  fluxes  for  each  surface  type  and  then  determines  grid- 
averaged  fluxes  based  on  the  fractional  coverage  of  each  surface  type  within  a 
grid  box. 

At  higher  resolution,  aggregation  of  fluxes  or  parameters  within  a  model 
grid  box  becomes  less  of  a  concern  than  the  advection  of  these  variables  from 
neighboring  grid  boxes,  particularly  when  surface  heterogeneities  are  at  least  as 
large  as  the  model  resolution  (personal  communication,  Wang  2010).  This  study 
is  intended  to  address  this  issue  by  examining  the  variation  of  surface  roughness 
with  wind  direction  near  a  coastline. 

C.  PURPOSE  AND  OVERVIEW  OF  THESIS 

In  this  thesis,  observational  data  from  a  dense  network  of  meteorological 
sensors  are  examined  to  explore  the  relationships  between  roughness  length  in 
the  lower  atmosphere  and  physical  properties  of  the  surface.  The  purpose  of  this 
study  is  to  demonstrate  a  local  dependence  of  the  downstream  turbulent 
response  to  varying  upstream  surface  features  at  different  length  scales  and  to 
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advocate  the  development  of  improved  algorithms  for  incorporating  this 
dependence  in  numerical  weather  and  climate  prediction  models.  The 
methodology  includes  analyzing  time  series  and  vertical  profiles  of  measured  and 
calculated  variables  in  order  to  evaluate  their  spatial  and  temporal  variations  and 
to  determine  their  relationships  to  surface  roughness  determined  from  satellite 
photographs. 

This  chapter  included  a  brief  introduction  of  terminology  and  concepts 
useful  in  this  study.  Chapter  II  presents  essential  concepts  of  thermodynamics  of 
the  lower  atmosphere  and  introduces  surface-layer  similarity  theory  and  flux- 
profile  relationships.  Chapter  III  describes  the  method  of  data  collection  and 
processing  and  the  meteorological  characteristics  of  the  region  for  the  cases 
presented.  Chapter  IV  details  the  results  and  analyses  from  the  selected  cases, 
and  Chapter  V  summarizes  the  results  and  concludes  the  thesis  with  a 
comparison  of  one  of  the  chosen  cases  to  numerical  model  output. 

D.  MILITARY  APPLICATIONS 

The  modern  military  operates  around  the  world  in  various  terrains  and 
land  surface  types,  often  in  regions  where  weather  observations  are  sparse.  The 
development  and  implementation  of  high-resolution  numerical  models  have  aided 
the  military  weather  community's  mission  to  enhance  operational  safety  while 
exploiting  the  weather  for  mission  success.  Communication  via  electromagnetic 
and  electro-optical  radiation,  low-level  aviation,  and  target  acquisition  and 
engagement,  are  all  examples  of  military  operations  that  are  directly  affected  by 
weather  conditions  in  the  near-surface  environment.  Turbulence  constantly 
alters  the  wind  field,  which  influences  the  distribution  of  visibility  restrictors  such 
as  low  clouds,  fog,  haze,  dust,  and  pollutants.  An  improved  understanding  of 
ABL  structure,  facilitated  by  studies  such  as  this  one,  will  lead  to  better 
techniques  for  simulating  the  effects  of  turbulent  processes  in  numerical  weather 
models  that  in  turn  lead  to  better  forecasts. 
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II.  SURFACE-LAYER  THERMODYNAMICS 


A.  TURBULENT  TRANSPORT 

Turbulence  describes  a  property  of  fluid  flows  characterized  by  strong 
spatial  gradients  in  velocity.  Turbulence  in  the  ABL  may  be  characterized  by  the 
presence  of  eddies,  complete  or  semi-complete  whorls  of  air  that  may  be 
transported  by  the  mean  wind  some  distance  downstream  from  their  point  of 
origin.  These  eddies  scale  on  a  spectrum  of  sizes  ranging  from  a  few  millimeters 
to  several  hundred  meters. 

Boundary  layer  meteorologists  typically  classify  turbulence  according  to 
the  process  that  generated  it.  Thermal  turbulence  occurs  when  parcels  of  air 
become  differentially  buoyant  by  the  addition  or  subtraction  of  thermal  energy, 
most  commonly  from  an  underlying  land  surface  heated  by  the  Sun  during  the 
daytime  or  cooled  at  night.  During  the  day,  these  parcels  or  “thermals” 
aggregate  near  the  surface,  forming  eddies  that  rise  semi-coherently  to  a  height 
where  their  densities  reach  equilibrium  with  that  of  the  surrounding  air,  at  which 
point  they  lose  upward  momentum  and  may  return  to  a  lower  level.  The  mean 
wind  may  also  transport  eddies  responding  to  this  thermally  induced  vertical 
motion  some  distance  from  their  point  of  genesis.  The  resultant  convective 
mixing  of  the  ABL  is  a  common  daytime  occurrence  over  land  in  temperate 
regions  throughout  the  world.  Mechanical  turbulence  in  the  ABL  is  associated 
with  swirling  air  currents  resulting  from  flow  around  surface  features  or  along  the 
interface  of  air  masses  with  distinct  densities.  In  general,  shear  increases 
(hence,  turbulence  is  stronger)  at  higher  wind  speeds  or  when  moving  air 
contacts  a  rougher  surface  (Stull  1988). 

Large  eddies  formed  by  either  thermal  or  mechanical  means  generate 
smaller  eddies  inertially  in  high-shear  regions  near  their  edges,  a  process  that 
effectively  transfers  turbulence  kinetic  energy  (TKE)  downscale.  Taken  together, 
these  variously  sized  eddies  comprise  an  energy  spectrum  for  turbulence.  TKE 


11 


contained  in  the  largest  eddies  “cascades”  through  the  inertial  subrange,  the 
domain  of  medium-sized  (about  102  m)  eddies.  Just  below  the  scale  of  the 
smallest  eddies  (about  10'3  m),  TKE  is  continually  dissipated  into  heat  by  fluid 
viscosity.  During  this  process,  TKE  is  not  a  conserved  atmospheric  property; 
rather,  conversion  of  TKE  to  internal  energy  is  constantly  occurring  at  the 
smallest  turbulent  scales  (McWilliams  2006).  Since  the  physical  properties  of  the 
underlying  surface  greatly  influence  the  production  and  evolution  of  both  thermal 
and  mechanical  turbulence  in  the  ABL,  models  attempting  to  approximate  the 
complexities  of  turbulent  transport  near  the  Earth’s  surface  must  incorporate 
accurate  information  about  these  properties. 

B.  SURFACE-LAYER  SCALING  PARAMETERS  AND  FLUXES 

A  frequently  used  scaling  parameter  for  turbulent  processes  in  the  surface 
layer  is  the  friction  velocity,  denoted  by  u,.  This  parameter  represents  the  effects 
of  surface  wind  stress  and  varies  with  both  wind  speed  and  surface  roughness 
(Kaimal  and  Finnigan  1994).  In  general,  u.t  is  given  by 

u *2  =  ylu'w'2  +  vV  .  (2) 

(Here,  u  ,  v' ,  and  w'  are  the  horizontal  and  vertical  components  of  air  velocity 
due  to  turbulence.)  For  simplicity,  ABL  meteorologists  typically  orient  their 
coordinate  axes  such  that  the  abscissa  corresponds  to  the  direction  the  wind 
stress  is  applied  (Stull  1988).  In  this  case,  one  horizontal  perturbation  dimension 
is  eliminated,  and  (2)  reduces  to 

-u4 2  =  u'w'  ,  (3) 

Where,  the  negative  sign  is  included  to  properly  orient  the  sign  of  momentum 
transfer.  Similarly,  a  scaling  parameter  for  surface-layer  virtual  potential 
temperature  (heat)  and  water  vapor  applications  is  introduced,  defined  by  the 
following: 
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Possibly  the  most  important  diagnosis  of  the  ABL  that  must  be  undertaken 
when  characterizing  turbulence  potential  is  that  of  its  static  stability. 
Meteorologists  usually  classify  a  layer’s  static  stability  somewhat  broadly  as 
stable,  unstable,  or  neutral.  (Here,  stability  is  qualified  with  the  word  “static”  to 
imply  that  this  atmospheric  condition  is  independent  of  air  motion.)  At  a  specific 
location,  if  air  with  a  particular  density  overlies  less  dense  air  (whether  by 
contrasting  temperature,  humidity,  or  both),  then  the  situation  is  unstable  and  the 
potential  for  spontaneous  convective  overturning  exists.  Instability  in  the 
boundary  layer  is  usually  a  result  of  heating  from  below,  although  many 
processes  tend  to  stabilize  the  boundary  layer,  including  radiative  emission  and 
molecular  diffusivity  (Stewart  1979).  Because  air  density  (at  constant  pressure) 
is  a  function  of  temperature  and  water  vapor  content,  absolute  virtual  potential 
temperature  is  used  for  calculations  involving  static  stability  in  this  study. 

When  quantifying  stability,  ABL  meteorologists  typically  employ  the 
concept  of  flux,  which  is  the  transport  of  a  variable  per  unit  area  per  unit  time 
(Stull  1988).  The  transport  of  thermal  energy,  referred  to  hereon  as  heat,  is  used 

to  demonstrate  the  concept  of  flux  used  by  fluid  dynamicists.  The  quantity  QH  is 
used  to  represent  the  transport  of  heat  by  a  fluid  medium  through  a  particular 
area  in  a  given  time.  Therefore,  in  SI  units,  QH  is  expressed  in  joules  per  square 
meter  per  second.  However,  atmospheric  heat  content  is  rarely  measured  in 
joules,  so  Qh  is  converted  to  the  more  convenient  kinematic  flux  form  QH  by 
dividing  the  flux  by  the  product  of  the  air  density  pair  and  the  specific  heat  of  air 
at  constant  pressure  Cp,  yielding 
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(6) 


Dimensional  analysis  of  the  result  reveals  the  principal  utility  of  this 
process;  namely,  that  the  quantity  QH  has  dimensions  of  absolute  temperature 

times  velocity,  both  of  which  are  routinely  measured  atmospheric  fields. 
Moreover,  except  in  extreme  cases,  the  density  of  the  lower  atmosphere  is 
roughly  constant  with  height,  typically  varying  by  a  factor  of  only  10%  from  the 
surface  to  the  top  of  the  ABL  (Stull  1988).  Using  this  assumption,  pajr  is  taken  to 

be  constant  rather  than  having  dependence  upon  pressure  and  temperature,  as 
in  the  equation  of  state  for  an  ideal  gas  (Bohren  and  Albrecht  1998).  To  simplify 
calculations  involving  observed  variables,  ABL  researchers  routinely  invoke 
similar  reasoning  in  the  literature;  and  hereinafter,  use  of  the  term  flux  implies 
kinematic  flux. 

C.  SURFACE-LAYER  FLUX-PROFILE  RELATIONSHIPS 

Though  our  mathematical  description  of  the  physics  is  imperfect,  many 
processes  associated  with  turbulence  in  the  ABL  show  remarkable  regularity 
over  a  wide  range  of  conditions  in  the  field,  suggesting  the  existence  of  functional 
relationships  between  the  variables  involved  (Stull  1988).  Historically,  ABL 
meteorologists  apply  dimensional  analysis  techniques  such  as  Buckingham  pi 
theory  to  investigate  these  relationships.  The  concept  of  Buckingham  pi  theory  is 
briefly  introduced  and  its  use  is  demonstrated  by  deriving  expressions  used  in 
calculations  for  this  study.  For  a  thorough  treatment  of  this  and  other  similarity 
techniques,  the  reader  is  referred  to  Kaimal  and  Finnigan  (1994). 

In  the  surface  layer,  strong  vertical  (and  horizontal)  time-dependent 
gradients  of  momentum,  mass,  and  energy  often  exist.  However,  the  associated 
vertical  fluxes  of  these  quantities  are  nearly  constant  with  height  (Kaimal  and 
Finnigan  1994).  The  objective  of  the  procedure  presented  below  is  to  determine 
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flux-profile  relationships,  empirical  associations  between  the  vertical  gradients  of 
these  measured  quantities  and  their  vertical  fluxes  due  to  turbulent  transport. 

First,  the  process  requires  identifying  variables  of  interest.  In  the  case  of 
turbulent  transfer  of  momentum,  the  variables  are:  z  ,  the  local  height  above  the 

dii 

surface;  — ,  the  vertical  gradient  of  the  mean  wind  u  ;  «*,  the  friction  velocity; 

dz 


turbulent  heat  flux.  Since  the  units  of  these  five  variables  comprise  three 
fundamental  dimensions — length,  time,  and  absolute  temperature — from  pi 
theory,  two  nondimensional  quantities  nx  and  n2  describe  the  transfer  of 

momentum  in  the  surface  layer.  Three  of  the  five  variables  are  selected  as  “key 
variables,”  that  together  contain  values  expressed  in  all  three  relevant 
dimensions.  The  remaining  two  variables  are  expressed  in  terms  of  the  key 
variables,  raised  to  the  unknown  real  exponents  a-f ,  as  follows: 


g_ 
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Examination  of  the  dimensions  of  the  variables  in  these  expressions  leads  to 
unique  solutions  for  a-f ,  yielding  the  dimensionless  quantities,  or  “groups,” 
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k  is  the  dimensionless  von  Karman  constant,  shown  by  analysis  of  field  data  to 
have  a  value  between  0.35  and  0.43  (Stull  1988)  (also  included  in  the  expression 
for  n2).  For  calculations  in  this  study,  the  most  common  approximation  for  the 
von  Karman  constant  is  adopted,  k  -  0.35. 

Note  that  the  sign  of  the  Monin-Obukhov  length  is  a  stability  indicator 
because  it  contains  the  turbulent  heat  flux  term  w'0'v ,  that  undergoes  a  sign 

change  when  static  stability  changes.  If  in  the  surface  layer  w'0'v  <  0,  thenZ  >  0 
and  the  heat  flux  is  positive  downward;  hence,  the  layer  is  stable.  Conversely,  if 
w'6'v  >  0,  heat  flux  is  positive  upward,  L  <  0,  and  the  layer  is  unstable.  If  there  is 

no  vertical  turbulent  heat  flux  (though  turbulence  can  be  present),  w'6'v  =  0,  Z 
approaches  (minus)  infinity,  and  the  layer  is  statically  neutral. 

From  (9)  and  (11),  nx  =  — .  According  to  pi  theory,  there  exists  a  functional 

Z 

dependence  ^  (though  not  necessarily  a  unique  one)  between  nx  and  n2_ 
expressed  by  the  nondimensional  gradient  relationship 

KZ  du  .  /  N.  ,  (  z\ 

X2=——  =  <PM^)  =  (t)m  —  >  (13) 

u,  dz  \L J 

with  the  subscript  m  used  to  indicate  momentum.  A  similar  procedure  yields  the 
gradient  relationships  for  the  mean  virtual  potential  temperature  profile 
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and  the  mean  specific  humidity  profile 


f  z  ^ 


kz  dq 
q ,  dz 


(15) 


where  6yt  and  qt  are  given  by  (4)  and  (5). 


Pioneering  research  by  Businger  et  al.  (1971)  and  Dyer  (1974)  involving 
analysis  of  surface-layer  observations  over  flat  terrain  in  varying  stability 
conditions  yielded  the  following  estimates  for  the  momentum  and  heat  flux-profile 
relationships  <f>m  and  <f>h q .  (The  flux-profile  relationship  for  specific  humidity  and 

other  scalar  fields  have  the  same  functional  form  as  that  for  heat.)  These  widely 
accepted  empirical  forms  are  used  to  analyze  the  data  in  this  study: 
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If  measurements  of  the  mean  wind  and  absolute  virtual  potential 
temperature  are  available  at  two  levels  z,  and  z2,  (1 3)— (1 5)  are  integrated  using 
these  limits  to  obtain 
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q*  = 
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(23) 


Although  the  integrals  in  (21)  -  (23)  are  expressible  in  closed  form,  they  must  be 
solved  iteratively  to  obtain  «*,  #„»,  or  q, ,  since  the  integrand  includes  Z  that 
involves  these  unknowns.  For  calculations  involving  (20)  -  (22)  undertaken  in 
this  study,  for  a  reasonable  initial  approximation  of  ut,  6vt ,  or  g*,  100  iterations 
converged  at  far  greater  precision  than  the  data  used  in  the  calculation. 

If  measurements  of  the  mean  wind  and  virtual  potential  temperature  are 
available  at  only  one  level  z,  (21)  -  (23)  are  supplemented  with  properties  from 
the  roughness  height  z0;  namely,  over  open  ocean  the  Charnock  relation  (1)  and 

the  assumptions  that  0sea  =^(z0),  qs(0sea)  =  q(zo)  and  in  general,  that  i7(z0)  =  0. 
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With  these  substitutions,  the  integrals  in  (21)  -  (23)  are  then  evaluated  from  z0  to 
z  to  obtain  an  expression  for  the  mean  wind  as  a  function  of  height,  from  which 
z0  can  be  determined. 


D.  BULK  AERODYNAMIC  PARAMETERIZATIONS 

In  practice,  most  numerical  models  incorporate  the  effects  of  surface 
roughness  by  parameterizations  based  on  the  bulk  aerodynamic  formulae.  The 
following  is  a  brief  introduction  of  the  concept,  adapted  from  Wang  (2009). 

The  drag  coefficient  C,  at  height  z  is  defined  as  follows: 


where  Uz  is  the  mean  wind  at  height  z  .  Researchers  working  within  the  surface 

layer  generally  adopt  a  standard  reference  height  (or  the  height  at  which  they 
measure  the  wind).  At  the  World  Meteorological  Organization’s  (WMO)  standard 
10-m  height  for  wind  observations,  C.  is  designated  CD.  In  (24),  this  becomes 


CD  = 


V^io  J 


Substituting  Equation  (3),  the  turbulent  momentum  flux  is  given  by 


u'w'  =  -CDUl02  =  -u2 

and  using  Equation  (21),  CD  is  given  by 


(25) 


(26) 


(27) 
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Similarly,  the  heat  exchange  coefficient  CH  (at  10  m)  is  given  by 
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III.  OBSERVATIONS  AND  SELECTED  CASES 


A.  SURFACE-LAYER  OBSERVATIONS 

The  National  Aeronautics  and  Space  Administration  (NASA)  Headquarters 
Weather  Support  Office  and  Cape  Canaveral  Air  Force  Station  (CCAFS)  jointly 
operate  the  Weather  Information  Network  Display  System  (WINDS),  part  of  an 
extensive  meteorological  observation  program  on  and  in  the  region  surrounding 
CCAFS  and  the  Kennedy  Space  Center  (KSC)  in  Florida.  In  continuous 
operation  since  1986,  the  WINDS  collects,  archives,  and  disseminates 
observations  measured  by  a  suite  of  over  200  wind,  temperature,  humidity,  and 
pressure  sensors  attached  to  46  instrumented  towers  distributed  across  a 
1,200-km2  region  (Figure  3).  The  primary  objective  of  the  WINDS  is  to  provide 
45th  Weather  Squadron,  NASA  Safety,  and  Range  Safety  personnel  with  a 
comprehensive,  real-time  description  of  weather  conditions  from  the  surface  to 
150  m  in  support  of  the  KSC  spaceflight  mission  (Computer  Sciences  Raytheon 
[CSR]  2006).  This  study  utilizes  observations  recorded  at  5-minute  intervals  from 
selected  WINDS  towers. 
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Figure  3.  WINDS  tower  distribution  on  and  near  CCAFS  on  the  Atlantic  Coast 
of  Florida.  After  CSR  (2006). 


WINDS  towers  at  CCAFS  are  classified  into  three  categories  based  on 
their  primary  operational  function:  launch  critical,  safety  critical,  and  forecast 
critical  towers.  The  four-digit  tower  identification  number  is  decoded  (except  for 
Towers  9001  and  9404)  as  follows:  the  first  pair  of  digits  is  the  tower’s  distance 
from  the  outer  coastline,  and  the  second  pair  of  digits  is  its  latitudinal  distance 
from  Port  Canaveral  (red  star  in  Figure  3).  Both  coded  distances  are  rounded  to 
the  nearest  integer  in  nautical  miles.  The  Appendix  contains  tables  of  tower 
locations  and  instrumentation  heights  by  tower  type. 


1.  Launch  Critical  Towers 

The  four  launch  critical  towers  are  located  nearest  the  active  launch 
complexes  and  the  Space  Shuttle  Landing  Facility.  Since  mission  planners  use 
observations  from  launch  critical  towers  to  directly  evaluate  launch  safety  criteria, 
their  accuracy  and  reliability  requirements  are  the  highest  in  the  WINDS.  Towers 
0002,  0006,  0110,  and  0313  all  contain  redundant  instrumentation  at  multiple 
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levels  and  are  equipped  with  battery  backup  to  ensure  data  availability  for  up  to 
24  hours  following  a  power  failure.  Horizontal  supports  for  wind  equipment  and 
temperature  sensors  are  aligned  northwest  to  southeast  on  Towers  0002,  0006, 
and  0110  and  northeast  to  southwest  on  Tower  0313  (Figure  4). 


Figure  4.  Photograph  of  the  upper  levels  of  Tower  0313,  a  150-m  launch 
critical  mast.  Note  the  dual  instrumentation  on  horizontal  supports 
at  each  level.  From  CSR  (2006). 

2.  Safety  Critical  Towers 

Fourteen  safety  critical  towers  are  primarily  located  near  facilities  where 
toxic  materials  such  as  vehicle  propellants  are  stored  (CSR  2006).  Safety  critical 
towers  support  wind  equipment  at  3.7  m  and  16.5  m  and  temperature  sensors  at 
1.8  m  and  16.5  m.  Since  relative  humidity  is  not  measured  at  these  locations, 
the  potential  temperature  0  is  used  in  calculations  involving  observations  from 
safety  critical  towers  instead  of  virtual  potential  temperature  6V .  Figure  5  is  a 
photograph  of  a  typical  safety  critical  tower. 
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Figure  5.  Photograph  of  a  typical  safety  critical  tower.  From  CSR  (2006). 

3.  Forecast  Critical  Towers 

Nineteen  forecast  critical  towers  are  primarily  distributed  in  the  region 
surrounding  CCAFS  (Figure  3).  Observations  from  these  towers  are  used  for 
general  weather  support  and  resource  protection  functions.  Like  launch  critical 
and  safety  critical  towers,  forecast  critical  towers  house  temperature  sensors  at 
1.8  m,  but  they  report  wind  data  from  one  level  only  (16.5  m)  (CSR  2006). 
Because  forecast  critical  towers  lack  multilevel  observations  of  wind  and 
temperature,  observations  recorded  at  these  locations  were  not  used  in  this 
study. 


4.  Data  Processing 

For  this  study,  observations  were  provided  with  the  following  precisions: 
air  temperature  and  dew  point,  0.1  °F;  surface  pressure  (Tower  0313  only), 
0.1  hPa;  wind  speed,  1  knot;  and  wind  direction,  1°.  Prior  to  using  these 
observations  to  derive  roughness  length,  friction  velocity,  and  heat  fluxes,  the 
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data  were  converted  to  corresponding  SI  units  and  subjected  to  quality  control 
algorithms  that  excluded  missing  and  unphysical  values.  Additionally,  since  the 
observation  interval  is  5  minutes,  all  variables  were  averaged  over  20  minutes  in 
order  to  more  accurately  represent  mean  conditions.  Due  to  the  relatively  coarse 
precision  of  wind  speed  measurements,  observations  taken  when  speed 
differences  between  the  two  levels  did  not  exceed  1  knot  (about  0.5  ms'1)  were 
excluded  from  calculations  of  z0 ,  u, ,  and  heat  fluxes. 

Figures  6  and  7  are  time  series  of  observed  conditions  and  calculated 
quantities  for  the  northeast  sensors  at  the  lowest  levels  (referred  to  in  this  study 
as  levels  1  and  2)  of  Tower  0313  on  2-5  May  2008.  Horizontal  axes  for  all  time 
series  presented  herein  are  labeled  with  the  day  of  the  year.  Local  midnight 
corresponds  to  the  day  number  tick  mark,  and  local  noon  coincides  with  an 
unlabeled  tick  mark.  Regardless  of  the  date,  time  is  given  in  Eastern  Standard 
Time  (EST)  without  correcting  for  Daylight  Savings  Time. 

In  Figures  6a  and  6b,  virtual  potential  temperature  (0v)  and  specific 
humidity  (q)  are  shown  for  the  1.8-m  (red)  and  16.5-m  (blue)  sensors.  Figures 
6c  and  6d  contain  time  series  of  the  calculated  sensible  heat  flux  (SHF)  and 
latent  heat  flux  (LHF).  For  the  flux  plots,  thermally  unstable  (red)  and  stable 
(blue)  conditions  are  determined  by  the  0v  difference  between  the  indicated 
levels.  Diurnally-cycling  stability  fluctuations  and  low-level  moisture  redistribution 
are  evident  in  the  0V  and  q  trends,  and  both  phenomena  appear  reasonable  in 

magnitude  and  duration.  Similarly,  fluxes  during  peak  afternoon  heating  fall 
within  normal  ranges  for  the  latitude  and  season  (personal  communication, 
Wang  2010). 
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Figure  6.  Temporal  variation  of  observed  (a)  absolute  virtual  potential 
temperature  (9V)  and  (b)  specific  humidity  (q)  in  g  kg'1  and 

calculated  (c)  sensible  and  (d)  latent  heat  fluxes  in  Wm'2  from  the 
northeast  sensors  at  temperature  levels  1  and  2  of  Tower  0313  on 
2-5  May  2008. 

Figures  7a  and  7b  contain  the  wind  speed  (WS)  and  wind  direction  (WD) 
observations  for  the  3.7-m  (red)  and  16.5-m  (blue)  sensors.  (Note  that  the  height 
of  the  lowest  wind  observation  differs  from  that  of  temperature  and  humidity  [cf. 
Appendix],  The  designation  "level  1"  is  used  herein  to  refer  to  both  the  lowest 
levels  for  wind  and  temperature.  Context  clarifies  whether  "level  1"  refers  to  the 
lowest  level  of  wind  or  temperature  measurement.)  On  days  122-124, 
progression  of  the  synoptic-scale  pattern  induced  a  gradual  veering  from 
southeasterly  to  southwesterly  flow,  with  higher  wind  speeds  observed  during  the 
daylight  hours  than  at  night.  Figures  7c  and  7d  are  the  roughness  length  (z0) 

and  friction  velocity  ( u *)  calculated  from  observations  at  the  indicated  levels  and 

shown  in  red  or  blue  based  on  the  same  convention  for  stability  as  in  Figures  6c 
and  6d. 
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Eastern  Standard  Time  (days) 

Figure  7.  Temporal  variation  of  observed  (a)  wind  speed  (WS)  in  ms'1  and  (b) 
wind  direction  (WD)  in  degrees  and  calculated  (c)  roughness  length 
(z0)  in  m  and  (d)  friction  velocity  (w«)  in  ms'1  from  the  northeast 
sensors  at  wind  levels  1  and  2  of  Tower  0313  on  2-5  May  2008. 

During  most  of  this  four-day  period,  z0  <  0.5  m,  and  w*  ranges  between 
0.2  and  0.6  ms'1.  However,  between  midnight  and  sunrise  of  day  123  and  for 
most  of  the  first  half  of  day  124,  unreasonably  large  z0  values  are  calculated. 
During  both  of  these  periods,  relatively  low  wind  speeds  (<  2  ms'1)  are  observed. 
In  addition,  large  temporal  variations  in  z0  are  noted  during  these  periods.  Similar 

trends  in  roughness  length  are  observed  throughout  the  year  during  periods  of 
low  wind  speeds,  often  occurring  overnight.  Calculated  roughness  lengths  during 
these  periods  are  considered  unphysical  for  two  reasons.  First,  the  wind  speed 
precision  (0.5  ms'1)  is  a  relatively  high  percentage  of  the  wind  speed.  The 
second  reason  is  that  Monin-Obukhov  similarity  theory  is  not  applicable  in  low 
wind  speed  conditions,  especially  for  nocturnal  stable  boundary  layers  in  which 
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turbulence  is  weak  and  intermittent.  For  this  reason,  in  the  remainder  of  the  time 
series  shown  herein,  results  for  z0are  not  plotted  if  the  wind  speed  at  either  level 
is  less  than  2  ms'1. 

Gaps  in  the  time  series  of  measured  fields  (seen  before  noon  on  day  122 
and  near  midnight  on  day  123)  correspond  to  periods  for  which  no  measurement 
is  available,  although  this  does  not  usually  indicate  a  sensor  has  failed.  Rather, 
the  problem  appears  to  occur  in  either  the  data  acquisition  or  storage  system  for 
the  particular  level  or  tower,  since  these  gaps  usually  occur  simultaneously  for 
both  wind  and  temperature  sensors,  often  at  multiple  levels.  Gaps  in  the  derived 
quantities  signify  at  least  one  of  the  following  conditions:  (1 )  lack  of  one  or  more 
measurements  required  for  the  calculation,  (2)  failure  to  exceed  the  0.5  ms'1 
speed  shear  criterion  between  the  indicated  levels,  (3)  wind  speed  at  either  level 
is  below  2  ms'1,  or  (4)  wind  speed  at  the  lower  level  exceeds  that  of  the  higher 
level. 

B.  CASES  SELECTED  FOR  THIS  STUDY 

1.  Surface  Morphology 

Located  on  the  East  Coast  of  the  Florida  Peninsula,  the  region 
surrounding  CCAFS  contains  diverse  types  of  surfaces  including  urban  areas, 
marshes,  temperate  forests,  barrier  islands,  sheltered  inland  waters,  and  the 
Atlantic  Ocean.  Depending  upon  the  strength  and  direction  of  the  surface-layer 
wind  field,  turbulent  eddies  generated  by  flow  over  these  heterogeneous  surfaces 
vary  in  dimension  and  impart  varying  gustiness  downstream. 

The  most  conspicuous  disparity  in  surface  roughness  at  CCAFS  exists 

between  land  and  ocean.  In  order  to  investigate  the  turbulent  response  that 

arises  from  this  heterogeneity,  periods  are  identified  during  which  near-surface 

winds  originated  over  land  (westerly)  and  ocean  (easterly).  These  opposite  flow 

regimes  should  be  selected  closely  in  time  to  minimize  the  effect  of  variations  in 

atmospheric  stability,  cloud  and  precipitation  forcing,  and  foliation.  Within  this 

framework,  this  study  focused  on  examining  the  turbulent  response  at  WINDS 
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towers  near  the  coast  during  a  diurnally  cycling  coastal  wind  circulation  event  in 
June  2008  and  the  passage  of  Tropical  Storm  Fay  in  August  2008.  Wind 
direction  changes  of  about  180°  occur  within  10-20  minutes  in  the  coastal  wind 
cases  and  in  less  than  12  hours  during  the  tropical  storm. 

2.  Coastal  Wind  Circulations 

Occurring  at  any  time  of  the  year,  but  most  prevalent  from  late  spring 
throughout  the  warm  season,  differential  heating  of  land  and  water  surfaces  near 
CCAFS  gives  rise  to  diurnally  cycling  coastal  winds,  commonly  referred  to  as  sea 
breezes  during  the  daytime  and  land  breezes  at  night.  The  leading  edge  of  the 
sea  breeze  defines  a  mesoscale  thermodynamic  discontinuity  called  the  sea- 
breeze  front  (SBF),  the  passage  of  which  is  usually  accompanied  by  a  slight 
decrease  in  temperature  and  a  sharp  increase  in  humidity.  During  periods  of 
relatively  weak  synoptic  forcing,  the  afternoon  SBF  may  travel  several  tens  of 
kilometers  inland  and  trigger  numerous  showers  and  thunderstorms  over  the 
Peninsula. 

Numerous  studies  of  coastal  wind  circulations  at  Cape  Canaveral  have 
been  conducted  in  situ  using  Doppler  radar  and  cloud  photogrammetry 
(Wakimoto  and  Atkins  1993);  observations  (Reed  1979);  aircraft  soundings  (Laird 
et  al.  1995);  and  with  high-resolution  mesoscale  modeling  (Rao  and  Fuelberg 
2000),  (Manobianco  et  al.  1996),  and  (Baker  et  al.  2000).  In  particular,  Reed 
(1979)  demonstrated  from  observations  at  Tower  0313  that  a  discernible  diurnal 
oscillation  between  onshore  and  offshore  wind  components  was  present  year- 
round,  with  the  largest  amplitude  occurring  in  May  and  the  smallest  in  January. 
Reed  also  found  the  largest  amplitudes  were  observed  by  the  highest 
anemometer  (then  151  m)  and  that  with  few  exceptions,  the  sea/land-breeze 
circulation  shifted  clockwise  (veered)  with  time,  roughly  completing  a  circuit  of  the 
compass  in  one  day. 
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In  Figures  6  and  7,  green  arrows  indicate  the  passage  of  two  SBFs  around 
noon  on  days  124  and  125  (4-5  May  2008).  SBF  passage  is  most  readily 
indicated  (in  this  region)  by  the  sudden  backing  of  the  wind  direction  from 
southwest  to  southeast  and  the  rapid  increase  in  q .  In  addition,  wind  speeds  at 
the  16.5-m  sensor  increase  from  about  2  ms'1  to  6  ms'1,  and  temperatures  at 
both  levels  gradually  fall  following  SBF  passage.  The  response  of  roughness 
length  to  coastal  wind  circulations  is  discussed  in  Chapter  IV. 

Figure  8  contains  vertical  profiles  of  0V,  q ,  wind  speed,  and  wind  direction 

observed  from  the  northeast  sensors  of  Tower  0313  at  4-hour  intervals  beginning 
at  0500  EST  on  4  May  2008.  Such  plots  are  useful  for  evaluating  the  evolution 
and  vertical  variation  of  these  measured  fields.  As  with  time  series  presented  in 
this  study,  observed  quantities  in  vertical  profiles  are  averaged  over  a  20-minute 
window  centered  on  the  indicated  time.  In  Figure  8a,  at  0500  EST  (about  1  hour 
prior  to  sunrise),  the  surface  layer  is  thermally  stable.  With  daytime  heating,  the 
lowest  layers  become  increasingly  unstable.  Around  1300  EST,  a  SBF  passed 
the  tower  location.  In  Figures  8c  and  8d,  the  wind  profiles  at  1300  EST  indicate 
that  the  SBF  had  passed  the  sensors  at  wind  levels  1-4,  but  had  not  yet  passed 
the  higher  levels  of  the  tower.  The  uniform  directional  profile  shown  before  and 
after  SBF  passage  was  typically  observed  on  days  affected  by  these  coastal 
wind  circulations.  In  Figure  8c,  logarithmic  wind  profiles  that  fit  to  observations 
from  the  lowest  two  levels  are  plotted  on  the  wind  speed  profiles  in  blue  dashed 
lines.  These  idealized  profiles  are  constructed  for  the  neutral  stability  condition 
using  the  methodology  presented  in  Chapter  II. 
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Tower  0313  (NE),  4  May  2008 


Figure  8.  Vertical  variation  of  observed  (a)  virtual  potential  temperature  (6V), 
(b)  specific  humidity  (q),  (c)  wind  speed  (WS),  and  (d)  wind 
direction  (WD)  from  the  northeast  sensors  of  Tower  0313  on  4  May 
2008. 


3.  Tropical  Storm  Fay 

The  easterly  wave  that  spawned  Tropical  Storm  Fay  moved  off  the  coast 
of  Africa  on  6  August  2008  and  rapidly  crossed  the  Atlantic  Ocean  before  slowing 
and  becoming  better  organized  over  the  Greater  Antilles.  On  15  August,  Fay 
was  upgraded  to  tropical  storm  status  while  over  land  on  the  island  of  Hispaniola. 
Under  favorable  upper-level  wind  conditions  and  despite  interacting  with 
mountainous  terrain,  Fay  continued  to  strengthen  slowly  and  made  two  landfalls 
in  Cuba  before  turning  northward  toward  the  Florida  Keys  (Stewart  and  Beven 
2009).  On  19-21  August,  Tropical  Storm  Fay  crossed  the  central  Florida 
Peninsula  before  making  its  third  Florida  landfall  near  Flagler  Beach  about  120 
km  north  of  CCAFS  (Figure  9).  While  near  CCAFS,  the  system’s  forward  speed 
slowed  to  3-4  knots,  resulting  in  16.5-m  sustained  winds  averaging  above  8  ms'1 
for  48  hours  at  most  WINDS  locations. 
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Figure  9.  Track  of  the  center  of  Tropical  Storm  Fay.  From  Stewart  and 
Beven  (2009). 


Tropical  Storm  Fay  proved  to  be  an  ideal  case  for  this  study  for  several 
reasons.  First,  the  track  of  Fay  caused  winds  to  shift  within  12  hours  from  a 
strong  onshore  to  a  strong  offshore  component.  Also,  Fay’s  prolonged  residence 
time  near  CCAFS  resulted  in  longer-fetch  surface-layer  winds  than  would  be 
induced  by  smaller-scale  flow  regimes.  Additionally,  fewer  wind  observations 
were  excluded  from  calculations  during  Fay  due  to  the  minimum  speed  shear 
criterion  of  0.5  ms'1.  This  is  a  result  of  both  the  relatively  low  ratio  of  this 
threshold  to  higher  wind  speeds  and  to  the  tendency  of  stronger  wind  profiles  to 
exhibit  significant  shear  in  the  surface  layer  over  land.  In  such  strong  wind 
profiles,  generation  of  ABL  turbulence  by  low-level  wind  shear  dominates  thermal 
turbulence  production,  and  surface-layer  stability  becomes  almost  thermally 
neutral.  Consequently,  uncertainties  associated  with  quantities  dependent  on 
thermal  stability,  such  as  those  in  the  Monin-Obukhov  similarity  functions  (16) 
and  (17),  are  minimized. 
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Figure  10  shows  time  series  of  observed  conditions  and  calculated 
variables  from  the  southeast  sensors  of  Tower  0006  during  the  passage  of 
Tropical  Storm  Fay  on  18-21  August  2008.  Variation  in  all  quantities  decreases 
on  days  232-233  due  to  clouds  and  precipitation  associated  with  the  tropical 
storm. 


Figure  10.  Same  as  Figure  6,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  temperature  levels  1 
and  2  of  Tower  0006  on  18-21  August  2008. 
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Tower  0006  (SE) 


Figure  1 1 .  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  1  and  2  of 
Tower  0006  on  18-21  August  2008. 


Figures  11a  and  11b  contain  time  series  of  the  wind  speed  and  direction 
observed  at  Tower  0006  during  Tropical  Storm  Fay.  While  the  center  of  Fay  was 
located  south  of  CCAFS  over  the  Florida  Peninsula  on  days  230  and  231,  flow 
was  generally  from  the  east.  As  Fay  approached  Cape  Canaveral,  wind  speeds 
increased,  reaching  a  maximum  of  about  16  ms'1  at  16.5  m  on  day  232.  The 
decrease  in  wind  speed  during  the  morning  hours  of  day  232  occurred  while  the 
center  of  Fay  was  near  CCAFS  (cf.  Figure  9).  During  about  a  12-hour  period,  the 
wind  direction  veered  to  the  west  and  speeds  increased  as  Fay  moved  north  of 
CCAFS.  The  response  of  roughness  length  during  Tropical  Storm  Fay  is 
discussed  in  Chapter  IV. 

Figure  12  contains  vertical  profiles  of  observed  conditions  recorded  by  the 
southeast  sensors  of  Tower  0006  at  4-hour  intervals  beginning  at  0400  EST  on 
20  August  2008.  In  Figure  12a,  the  near-neutral  thermal  stability  condition  is 
evidenced  by  the  nearly  vertical  0y  profiles.  Figures  12c  and  1 2d  depict  the  wind 


34 


profiles  before  and  after  the  center  of  Fay  crossed  the  latitude  of  the  tower,  with 
southeasterly  flow  at  0400  EST  and  0800  EST  and  southwesterly  flow  thereafter. 


Figure  12.  Same  as  Figure  8,  except  observed  conditions  are  from  the 
southeast  sensors  of  Tower  0006  on  20  August  2008. 
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IV.  RESULTS  AND  ANALYSES 


A.  ROUGHNESS  LENGTH  VARIATION  IN  RESPONSE  TO  COASTAL 

WIND  CIRCULATIONS 

On  1-4  June  2008,  a  series  of  coastal  wind  circulations  occurred  at 
CCAFS,  and  these  were  recorded  at  WINDS  locations  near  the  coast.  Daily 
mesoscale  sea  breezes  are  a  common  occurrence  in  the  region  throughout  the 
warm  season,  but  this  period  was  particularly  selected  for  its  lack  of  significant 
convective  activity  in  the  vicinity  of  CCAFS.  Although  Monin-Obukhov  surface- 
layer  similarity  can  be  applied  to  convective  winds,  precipitation-induced  outflows 
often  disrupt  the  repetitive  signals  in  observations  of  the  diurnal  coastal  wind 
cycle,  and  it  is  this  regularity  this  study  sought  to  investigate. 

1.  Tower  0006 

Tower  0006  is  a  launch  critical  tower  located  approximately  600  m  from 
the  shoreline  (Figures  3  and  13).  At  higher  magnification  (Figure  14),  the 
landscape  within  a  200-m  radius  of  Tower  0006  consists  of  mainly  shrubs  and 
small  trees,  the  largest  of  which  (using  a  qualitative  shadow  analysis  technique) 
lie  in  sectors  040°-120°,  150°-190°,  and  330°-350°.  (Herein,  all  sectors  and 
radials  are  presented  with  reference  to  the  tower  and  are  estimated  to  the 
nearest  10°  using  satellite  imagery.)  Within  a  50-m  radius  of  Tower  0006,  an  "X"- 
shaped  area  of  cleared  ground  contains  the  attachments  for  the  tower's  guy 
wires.  Variations  of  this  clearing  are  present  at  each  launch  critical  tower. 
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100  m  m 

Figure  13.  Satellite  photograph  of  Tower  0006  and  surrounding  area.  Red  "X" 
marks  the  tower  location.  In  all  satellite  photographs  presented 
herein,  true  north,  or  360°,  is  oriented  toward  the  top  of  the  image. 
Image  ©2010  Google. 


Figure  14.  Satellite  photograph  of  Tower  0006  and  immediate  surroundings. 
Red  "X"  marks  the  tower  location.  Image  ©2010  Google. 
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Figures  15  and  16  are  time  series  showing  observations  and  calculated 
quantities  from  the  southeast  sensors  at  levels  1  and  2  of  Tower  0006  on  1-4 
June  2008.  In  Figure  15a,  the  diurnal  temperature  variation  is  observed  at  both 
levels.  The  magnitude  of  the  temperature  change  is  greater  for  the  lower  level. 
In  Figure  15b,  specific  humidity  decreases  due  to  vertical  mixing  during  the 
morning  hours  then  rises  rapidly  in  the  marine  air  following  each  SBF.  In  Figures 
15c  and  15d,  diurnal  variation  in  surface-layer  thermal  stability  and  in  the 
calculated  SHF  and  LHF  values  are  also  apparent.  Maximum  values  for  SHF 
and  LHF  during  the  daytime  reach  about  300  and  200  Wm'2,  respectively. 

In  Figures  16a  and  16b,  wind  speed  and  direction  indicate  the  passage  of 
four  SBFs  before  noon  each  day  (indicated  throughout  by  green  arrows). 
Following  each  SBF,  wind  speeds  increase  and  the  flow  backs  rapidly  from  the 
southwest  to  the  southeast,  followed  by  a  gradual  veering  to  the  southwest 
during  the  remainder  of  the  day  and  overnight.  Maximum  wind  speeds  of  about  5 
ms'1  at  both  levels  occur  following  the  SBF  each  afternoon,  and  the  weakest 
winds  are  observed  around  midnight.  The  calculated  friction  velocity  ut  ,  shown 

in  Figure  16d,  ranges  between  0.1  and  0.5  ms'1,  a  typical  range  for  surface 
layers  over  land  at  these  wind  speeds.  Figure  16c  contains  the  calculated 
roughness  length  z0,  which  is  mostly  less  than  0.25  m.  As  stated  in  Chapter  III, 

no  z0  value  is  plotted  when  the  wind  speed  at  either  level  is  less  than  2  ms"1. 

On  all  four  days,  the  lowest  z0  values  for  levels  1  and  2  were  recorded 

when  the  wind  direction  was  from  sector  120°-150°,  which  is  seen  in  Figure  14  to 
roughly  correspond  to  an  opening  in  the  trees  connecting  ground  in  the 
immediate  vicinity  of  Tower  0006  to  close-cut  grass  and  roads  near  buildings 
located  200  m  southeast  of  the  tower.  At  approximately  this  distance  in  any 
other  direction  from  the  tower,  more  numerous  low  shrubs  and  trees  are  present, 
and  when  the  flow  is  from  these  sectors,  higher  z0  values  are  calculated  for 
levels  1  and  2. 
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Eastern  Standard  Time  (days) 


Figure  15.  Same  as  Figure  6,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  temperature  levels  1 
and  2  of  Tower  0006  on  1-4  June  2008. 


Eastern  Standard  Time  (days) 


Figure  16.  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  1  and  2  of 
Tower  0006  on  1-4  June  2008. 


40 


Figure  17  is  similar  to  Figure  16,  except  quantities  are  from  wind  levels  2 
and  3  of  Tower  0006.  Since  temperature  and  humidity  measurements  are  not 
recorded  at  wind  level  3,  those  from  levels  1  and  2  (Figure  15)  are  used  to  obtain 
z0  and  ut  at  levels  2  and  3.  Also,  due  to  the  lack  of  temperature  sensors  at  level 

3,  heat  fluxes  calculated  using  wind  levels  2  and  3  and  temperature  levels  1  and 
2  at  are  very  similar  to  those  in  Figures  13c  and  13d  and  are  not  shown  for  this 
case.  Although  wind  speed  at  level  3  also  decreases  overnight,  winds  are 
generally  stronger  at  this  level,  most  notably  following  each  SBF. 


Eastern  Standard  Time  (days) 

Figure  17.  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  2  and  3  of 
Tower  0006  on  1-4  June  2008. 

Consistent  with  results  from  wind  levels  1  and  2,  at  wind  levels  2  and  3  the 
lowest  z0  values  occur  when  the  wind  direction  is  from  sector  120°-140°  on  days 

152-154.  Again,  this  response  appears  to  be  related  to  the  smoother  landscape 
in  this  sector  (cf.  Figure  13),  though  the  distance  to  the  surface  inducing  the  low 
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z0  values  in  this  layer  is  not  apparent.  With  the  exception  of  a  few  scattered 

trees  and  buildings,  this  sector  contains  a  relatively  smooth  landscape  up  to  500 
m  from  the  tower. 

Figure  17  also  contains  fewer  calculated  values  of  z0  and  u,  than  were 

possible  at  the  lowest  levels.  This  is  because  the  shear  between  levels  2  and  3 
during  this  period  more  frequently  failed  to  meet  the  minimum  criterion  of  0.5  ms' 
\  For  this  reason,  wind  levels  3  and  4  of  Tower  0006  (likewise  for  Tower  0110  in 
the  following  section),  that  contain  even  more  such  exclusions,  are  not  shown 
here  for  the  coastal  wind  cases. 

Figure  18  clearly  demonstrates  the  relationship  between  z0  and  wind 
direction  in  varying  stability  conditions  at  the  two  lowest  wind  "layers"  of  Tower 
0006.  It  contains  polar  plots  of  z0  averaged  over  15°  sectors  for  the  entire  year 

2008  for  the  southeast  sensors  at  levels  1  and  2  (Figure  18a)  and  levels  2  and  3 
(Figure  18b).  The  tower  location  is  the  pole  for  these  plots,  and  wind  direction  is 
the  azimuthal  coordinate.  Line  segments  connecting  the  average  z0  value 

(plotted  on  the  center  radials  of  each  15°  sector)  form  the  blue  outline  in  each 
plot.  Like  the  time  series  presented  in  this  study,  for  all  polar  plots  presented 
herein,  in  addition  to  averaging  all  variables  over  20  minutes,  z0  was  not 
calculated  for  wind  speeds  less  than  2  ms"1. 
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Figure  18.  Azimuthal  variation  of  calculated  roughness  length  z0  averaged 

over  each  15°  sector  from  the  southeast  sensors  at  (a)  wind  levels 
1  and  2  and  (b)  wind  levels  2  and  3  of  Tower  0006  for  1  January  to 
31  December  2008. 


These  "z0  roses"  compare  reasonably  well  to  the  satellite  images  of 
Tower  0006  (Figures  13  and  14).  In  Figure  18a,  the  "X"  pattern  is  discernible, 
though  the  branches  of  the  "X"  in  the  polar  plot  correspond  to  larger  z0  values 
associated  with  areas  of  trees  and  shrubs  between  the  cleared  areas  visible  on 
the  satellite  photograph.  In  Figure  18b,  sectors  with  larger  z0  correspond  to  wind 
directions  from  forested  areas  situated  100-300  m  from  the  tower.  At  the  same 
distance  to  the  west,  the  landscape  appears  to  be  smoother,  and  z0  values 
corresponding  to  westerly  winds  are  lower. 

2.  Tower  0110 

Tower  0110,  a  launch  critical  tower  housing  instrumentation  at  the  same 
levels  as  Tower  0006,  is  located  approximately  1.7  km  from  the  shoreline 
(Figures  3  and  37)  and  7  km  NNW  of  Tower  0006.  The  eastern  shore  of  the 
Banana  River  (actually  a  large  brackish  lagoon)  is  about  50  m  west  of  the  tower. 
At  higher  magnification  (Figure  20),  trees  are  visible  between  the  tower  and  the 
Banana  River,  and  trees  are  also  present  in  sectors  160°-200°  along  the 
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riverbank  and  across  Titan  III  Road  to  the  east.  Shadow  analysis  suggests  the 
largest  trees  in  the  immediate  vicinity  are  those  located  25-50  m  west-southwest 
of  the  tower  along  the  riverbank.  The  smoothest  sector  in  the  immediate  vicinity 
of  the  tower  is  that  containing  the  tower's  access  road  along  radial  100°  and  the 
southeast  branch  of  the  "X."  The  surface  in  this  sector  appears  relatively  smooth 
to  a  distance  of  about  75  m. 
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Figure  19.  Satellite  photograph  of  Tower  0110  and  surrounding  area.  Red  "X" 
marks  the  tower  location.  Image  ©2010  Google. 
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Figure  20.  Satellite  photograph  of  Tower  0110  and  immediate  surroundings. 

Red  "X"  marks  the  tower  location.  Image  ©2010  Google. 

Figures  21  and  22  are  time  series  showing  observations  and  calculated 
quantities  from  the  southeast  sensors  at  levels  1  and  2  of  Tower  0110  on  1-4 
June  2008.  As  observed  at  Tower  0006,  SBFs  occur  before  noon  each  day  with 
similar  responses  in  observed  variables.  In  Figures  22a  and  22b,  secondary 
circulations  (indicated  by  purple  arrows),  apparently  influenced  by  the  Banana 
River  (Laird  et  al.  1995),  are  observed  around  sunrise  each  day  as  about  a  45° 
veering  in  wind  direction  and  a  slight  increase  in  q.  Wind  speeds  decrease  at 
night  as  with  Tower  0006,  though  winds  at  Tower  0110  remain  stronger  on 
average  throughout  the  period,  possibly  enhanced  by  divergence  over  the 
Banana  River.  In  Figure  22c,  z0  values  less  than  0.25  m  are  calculated  in  the 
afternoons  following  the  SBFs,  during  which  periods  the  wind  direction  is  mostly 
from  sector  1 10°-140°.  As  at  Tower  0006,  the  z0  response  at  this  layer  appears 
to  be  due  to  the  smooth  surface  near  the  tower  in  that  direction. 
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Roughness  lengths  are  greater  (about  0.3-0. 5  m)  during  two  periods  of 
southerly  wind  at  speeds  of  3-5  ms'1,  one  just  before  midnight  on  day  153  and 
one  between  midnight  and  sunrise  on  day  155.  In  this  sector,  tree  density  is 
relatively  high,  and  the  larger  z0  values  appear  to  be  a  direct  result  of  flow  over 
this  denser  forest  canopy. 


Eastern  Standard  Time  (days) 


Figure  21.  Same  as  Figure  6,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  temperature  levels  1 
and  2  of  Tower  01 10  on  1-4  June  2008. 
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Eastern  Standard  Time  (days) 

Figure  22.  Same  as  Figure  7,  except  observed  conditions  and  calculated 


quantities  are  from  the  southeast  sensors  at  wind  levels  1  and  2  of 
Tower  01 10  on  1-4  June  2008. 


Figure  23  is  similar  to  Figure  22,  but  for  wind  levels  2  and  3  of  Tower 
01 10.  In  this  layer,  as  the  flow  veers  from  about  120°  to  170°  following  each  SBF 
passage,  z0  values  increase  from  about  0.5  m  to  1.5  m.  Two  very  brief  periods 

during  which  the  wind  speed  exceeded  5  ms'1  under  southwesterly  flow  are 
observed  just  after  sunrise  on  days  153  and  155  following  the  passage  of  the 
secondary  circulations  associated  with  the  Banana  River  (purple  arrows).  During 
these  periods,  z0  values  fall  to  less  than  0.2  m.  These  low  roughness  lengths 

are  an  apparent  response  to  flow  over  the  smooth  surface  of  the  Banana  River. 
Due  to  the  intervening  trees  west  of  the  tower  and  their  role  in  the  formation  of  an 
IBL,  the  Banana  River  was  not  "seen"  in  calculations  of  z0  at  levels  1  and  2  from 
any  wind  sector. 
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(a)  Tower  0110  (SE) 


Figure  23.  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  2  and  3  of 
Tower  01 10  on  1-4  June  2008. 


Figure  24  contains  "z0  roses"  for  Tower  0110  for  the  lowest  two  layers  of 
Tower  0110  for  the  entire  year  2008.  At  the  lowest  layer  (Figure  24a),  the 
familiar  branches  of  the  "X"  are  present,  and  a  lobe  of  z0  values  exceeding  2  m 

is  located  in  sector  250°-280°.  As  seen  in  Figure  20,  this  sector  contains  trees 
located  between  the  tower  and  the  Banana  River.  Between  the  16.5-m  and 
49.4-m  wind  levels  (Figure  24b),  the  IBL  formed  by  these  trees  is  not  as 
prevalent,  and  flow  from  this  sector  corresponds  to  lower  z0  values. 
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Figure  24.  Azimuthal  variation  of  calculated  roughness  length  z0  averaged 


over  each  15°  sector  from  the  southeast  sensors  at  (a)  wind  levels 
1  and  2  and  (b)  wind  levels  2  and  3  of  Tower  01 1 0  for  1  January  to 
31  December  2008. 

B.  ROUGHNESS  LENGTH  VARIATION  IN  RESPONSE  TO  TROPICAL 

STORM  FAY 

As  discussed  in  Chapter  III,  the  passage  of  Tropical  Storm  Fay  near 
CCAFS  on  20-21  August  2008  provided  an  ideal  set  of  observations  for 
consideration  in  this  study.  Launch  critical  towers  0002  and  0110  were  selected 
for  analysis  of  the  roughness  length  response  during  Tropical  Storm  Fay. 

1.  Tower  0002 

Tower  0002  is  located  approximately  700  m  from  the  coastline  (Figure  25). 
At  higher  magnification  (Figure  26),  the  surface  within  a  150-m  radius  of  the 
tower  is  fairly  uniform,  consisting  mainly  of  shrubs  and  small  trees  except  for  the 
familiar  "X"  pattern  and  an  effectively  elongated  southeast  branch  of  the  "X" 
aligned  with  the  tower's  access  road  along  radial  140°.  Qualitative  shadow 
analysis  reveals  the  largest  trees  in  the  area  around  Tower  0002  lie  about  200  m 
distant  in  sector  250°-330°.  At  a  radius  of  about  300-500  m,  the  surface  is 
somewhat  smoother,  both  along  the  immediate  shoreline  to  the  southeast  and  on 
cleared  land  containing  scattered  buildings  and  structures  to  the  west. 
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Figure  25.  Satellite  photograph  of  Tower  0002  and  surrounding  area.  Red  "X" 
marks  the  tower  location.  Image  ©2010  Google. 


80m 


Figure  26.  Satellite  photograph  of  Tower  0002  and  immediate  surroundings. 
Red  "X"  marks  the  tower  location.  Image  ©2010  Google. 

Figures  27  and  28  are  time  series  showing  measurements  and  calculated 
quantities  from  the  southeast  sensors  at  levels  1  and  2  of  Tower  0002  on  18-21 
August  2008.  In  Figure  27,  diurnal  variation  of  the  surface-layer  temperature  and 
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heat  fluxes  are  clearly  seen  on  day  230  but  become  obscured  on  the  remaining 
days  because  of  cloud  cover  and  precipitation  associated  with  Fay.  On  days 
231-234,  the  surface  layer  is  slightly  unstable,  SHF  averages  50  Wm'2,  and  LHF 
averages  100  Wm'2.  Near  the  end  of  the  period,  the  level-1  temperature  and 
humidity  sensors  failed  to  record  usable  data,  and  as  a  result,  calculations 
following  the  failure  should  be  disregarded. 


Eastern  Standard  Time  (days) 

Figure  27.  Same  as  Figure  6,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  temperature  levels  1 
and  2  of  Tower  0002  on  18-21  August  2008. 

In  Figure  28,  as  the  center  of  Tropical  Storm  Fay  approached  CCAFS  on 
days  230  and  231,  the  wind  direction  remained  fairly  consistent  from  the  east  and 
southeast  and  gradually  increased  in  speed.  With  appreciable  shear  between 
wind  levels  1  and  2,  u„  also  increased  with  wind  speed  as  expected.  On  day 

231,  with  a  steady  east  wind,  z0  values  averaged  0.5  m,  a  reasonable  magnitude 

given  the  numerous  but  separate  small  trees  in  that  sector  within  100-200  m  of 
the  tower  (Figure  26).  From  midnight  until  noon  on  day  232,  the  wind  direction 
veered  at  a  fairly  constant  rate  from  easterly  to  westerly  as  the  center  of  Fay 
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crossed  the  latitude  of  CCAFS.  During  this  12-hour  period,  with  wind  speeds 
exceeding  10  ms'1,  z0  values  sharply  decreased  to  about  0.1  m  when  the  wind 
direction  was  from  sector  120°-150°,  which  corresponds  closely  to  the  sector 
containing  the  tower  access  road.  By  about  sunrise  (0600  EST),  the  wind 
direction  had  shifted  more  to  the  south  and  z0  increased,  which  is  reasonable 

given  the  trees  present  in  that  sector.  Around  noon,  a  brief  period  of  low  z0 

values  is  identifiable  when  the  flow  is  from  about  radial  220°,  which  contains 
cleared  ground  along  the  southwest  branch  of  the  "X".  After  backing  about  60°, 
the  wind  remains  steady  from  about  220°  on  day  233  and  calculated  z0  values 

are  slightly  lower  as  a  result  of  the  relative  smoothness  of  this  sector  within  50  m 
of  the  tower. 
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Figure  28.  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  1  and  2  of 
Tower  0002  on  18-21  August  2008. 


Figure  29  shows  measurements  from  the  southeast  sensors  at  wind  levels 
2  and  3  of  Tower  0002  during  Tropical  Storm  Fay.  At  these  levels,  z0  (calculated 

using  thermal  data  from  levels  1  and  2)  averages  above  1  m  until  just  before 
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noon  on  day  232,  when  values  sharply  decrease.  During  this  brief  period  of  low 
z0  values  (about  0.2  m),  the  flow  is  from  sector  220°-240°.  In  addition  to  the  "X" 
branch  in  this  sector,  cleared  land  southwest  of  the  tower  is  closest  (about  220 
m)  in  this  direction  and  may  be  contributing  to  the  low  z0 .  If  this  is  true,  then  this 

observation  also  demonstrates  the  dependence  of  z0  upon  wind  speed  since 
lower  speeds  in  the  coastal  wind  circulation  cases  corresponded  to  z0  values 
indicative  of  shorter  distances  "seen"  by  wind  levels  2  and  3. 


Figure  29.  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  2  and  3  of 
Tower  0002  on  18-21  August  2008. 

Figures  30  and  31  are  from  the  southeast  sensors  at  temperature  levels  2 

and  3  and  wind  levels  4  and  5  of  Tower  0002  during  Tropical  Storm  Fay.  The 

large  gaps  in  the  calculated  values  on  day  230  are  due  to  insufficient  wind  shear. 

In  Figure  30,  specific  humidity  at  level  3  is  suspect,  since  the  trace  of  q  at  both 

levels  is  essentially  parallel  and  the  16.5-m  sensor  averages  2  g  kg'1  higher  than 

the  62.2-m  sensor.  Such  a  large  difference  in  specific  humidity  in  less  than  50  m 

seems  unrealistic  given  the  uniformity  of  the  wind  field.  On  day  233,  temperature 
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and  humidity  fails  to  record  at  level  3,  and  calculated  results  thereafter  are  not 
reliable.  Despite  this  questionable  humidity  data  on  days  230-233,  confidence  in 
z0  remains  high  since  the  effect  of  variations  in  stability  is  minimized  when  wind 
speed  is  stronger  in  near-neutral  conditions,  as  noted  in  Chapter  III. 

In  an  attempt  to  both  validate  the  hypothesis  of  this  study  and  quantify  the 
degree  to  which  the  distance  "seen"  by  various  tower  layers  varies  with  altitude 
and  wind  speed,  it  is  tempting  to  search  for  evidence  that  Tower  0002  can  "see" 
the  ocean  from  observations  at  its  highest  levels.  In  Figure  31 ,  from  noon  on  day 
231  to  midnight  on  day  232,  with  a  steady  easterly  wind,  z0  values  are  between  1 
m  and  3  m.  In  this  direction,  as  seen  in  Figure  32,  the  distance  to  the  ocean  is 
approximately  1.6  km.  Once  the  wind  direction  begins  to  veer  on  day  232,  z0 

decreases  to  less  than  0.5  m,  remaining  at  about  this  value,  while  the  flow  is  from 
sector  1 10°-180°.  Such  low  roughness  lengths  in  this  sector  could  be  caused  by 
either  the  ocean  (presumably  roughened  due  to  the  tropical  storm)  or  the  sparse 
vegetation  along  the  immediate  shoreline,  since  the  distance  between  the  tower 
and  the  ocean  is  lowest  in  this  sector. 
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Figure  30.  Same  as  Figure  6,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  temperature  levels  2 
and  3  of  Tower  0002  on  18-21  August  2008. 
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Figure  31.  Same  as  Figure  7,  except  observed  conditions  and 
calculated  quantities  are  from  the  southeast  sensors  at  wind  levels 
4  and  5  of  Tower  0002  on  1 8-21  August  2008. 


1km 


Figure  32.  Satellite  photograph  showing  the  proximity  of  Tower  0002  to  the 
Atlantic  Ocean.  The  tower  is  less  than  1  km  from  the  ocean  in 
sector  120°-180°.  Red  "X"  marks  the  tower  location.  Image  ©2010 
Google. 
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Figure  33  contains  polar  plots  of  z0  for  the  lowest  two  layers  of  Tower 

0002  for  the  entire  year  2008.  In  Figure  33a,  the  lowest  layer  is  again  dominated 
by  the  "X"  pattern,  which  is  more  uniform  here  than  for  Towers  0006  and  0110. 
In  Figure  33b,  z0  maxima  in  sector  260°-280°  and  along  radial  320°  correspond 


to  trees  about  200  m  from  the  tower  as  previously  noted.  The  decrease  in  z0 


observed  between  these  sectors  corresponds  to  a  gap  in  the  trees  visible  in  the 
upper  left  of  Figure  26,  which  opens  to  cleared  land  beyond.  Another  relative 
minimum  occurs  in  sector  220°-250°,  which  contains  cleared  land  200-300  m 
from  the  tower. 
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Figure  33.  Azimuthal  variation  of  calculated  roughness  length  z0  averaged 

over  each  15°  sector  from  the  southeast  sensors  at  (a)  wind  levels 
1  and  2  and  (b)  wind  levels  2  and  3  of  Tower  0002  for  1  January  to 
31  December  2008. 


2.  Tower  0110 

In  addition  to  the  coastal  wind  circulation  case,  observations  from  Tower 
0110  were  also  selected  to  evaluate  the  z0  response  during  Tropical  Storm  Fay. 

Satellite  photographs  of  the  tower  are  given  in  Figures  19  and  20.  For  this  case, 
the  extent  to  which  the  Banana  River  can  be  "seen"  by  the  various  levels  of 
Tower  0110  and  whether  this  effect  differs  from  the  coastal  wind  case  are 
investigated. 
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Figures  34  and  35  are  time  series  showing  measurements  and  calculated 
quantities  from  the  southeast  sensors  at  levels  1  and  2  of  Tower  0110  on  18-21 
August  2008.  In  Figure  34,  observations  of  temperature  and  humidity  and  the 
calculated  heat  fluxes  at  this  location  are  similar  to  that  of  Tower  0002,  located 
16  km  to  the  southeast.  Just  after  midnight  on  day  232,  the  level  2  humidity 
sensor  fails  to  record  data  but  appears  to  recover  on  day  233.  Again,  this  outage 
does  not  significantly  affect  the  calculations  for  z0  due  to  the  near-neutral 
conditions  during  this  event. 


Eastern  Standard  Time  (days) 

Figure  34.  Same  as  Figure  6,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  temperature  levels  1 
and  2  of  Tower  01 10  on  18-21  August  2008. 

In  Figure  35,  the  wind  direction  was  from  the  east  and  southeast,  and  the 
wind  speed  steadily  increased  until  the  morning  of  day  232.  Just  before  noon,  as 
the  center  of  Fay  was  near  the  tower,  winds  calmed  and  abruptly  veered  to  the 
west  in  about  1  hour,  a  notable  difference  from  the  12-hour  wind  shift  observed  at 
Tower  0002.  On  days  231  and  232,  lower  z0  values  are  calculated  while  the 
wind  direction  was  from  the  east  and  southeast,  a  similar  trend  to  the  coastal 
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wind  case.  After  the  wind  shift,  z0  rapidly  increases  to  over  2  m,  which  is 
reasonable  given  the  wind  speed  and  the  proximity  of  the  trees  west  of  the  tower. 
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Figure  35.  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  1  and  2  of 
Tower  01 10  on  18-21  August  2008. 


Figure  36  shows  observed  and  calculated  values  from  wind  levels  2  and  3 
of  Tower  0110  during  Tropical  Storm  Fay.  At  these  levels,  z0  (calculated  using 
thermal  data  from  levels  1  and  2)  rapidly  decreases  at  the  time  of  the  wind  shift. 
Corresponding  to  flow  from  the  east  and  southeast,  z0  values  averaged  above  1 

m,  while  after  the  shift,  z0  values  fall  to  0.05  m,  indicative  of  flow  over  the  very 
smooth  surface  of  the  Banana  River.  Near  the  end  of  the  period,  the  wind 
direction  has  backed  to  about  180°,  and  z0  values  increase  as  a  result  of  the 
fetch  over  forested  land.  Comparing  Figures  23  and  36,  it  is  possible  to  discern  a 
similar  effect  in  the  z0  response  at  wind  levels  2  and  3  in  both  the  coastal  wind 

and  Fay  cases,  although  fewer  observations  of  wind  from  the  direction  of  the 
River  were  recorded  for  the  coastal  wind  case. 
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Figure  36.  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  2  and  3  of 
Tower  01 10  on  18-21  August  2008. 


Figure  37.  Satellite  photograph  showing  the  proximity  of  Tower  0110  to  the 
Banana  River  and  Atlantic  Ocean.  Red  "X"  marks  the  tower 
location.  Image  ©2010  Google. 
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In  Figure  37,  the  west  to  southwest  wind  sector  following  the  shift  has  a 
1-2  km  fetch  across  the  Banana  River,  and  the  very  low  z0  values  in  Figure  35 

confirm  that  wind  levels  2  and  3  can  "see"  this  fetch,  while  levels  1  and  2  sense 
the  IBL  induced  by  the  trees  along  the  riverbank. 

Figures  38  and  39  are  time  series  of  observed  and  calculated  values  from 
temperature  levels  2  and  3  and  wind  levels  3  and  4  of  Tower  0110  during 
Tropical  Storm  Fay.  As  observed  at  the  upper  levels  of  Tower  0002,  lower  shear 
at  these  levels  results  in  more  gaps  in  the  calculated  quantities  on  both  figures. 
Although  relatively  few  data  points  for  z0  are  present  on  day  233,  they  are  of 

similar  magnitude  to  those  observed  at  wind  levels  2  and  3  and  support  the 
conclusion  that  these  low  roughness  lengths  result  from  flow  across  the  Banana 
River. 

Before  sunrise  on  day  232  and  near  the  time  of  the  maximum  observed 
wind  speed  of  25  ms'1,  a  second  period  of  low  z0  values  correspond  to  a  wind 

direction  of  about  120°.  Along  this  radial,  areas  of  trees  and  the  northern  edge  of 
cleared  land  around  the  Mars  Observer  Launch  Facility  (Figure  37,  lower  right)  lie 
between  the  tower  and  the  ocean,  which  is  about  2  km  from  the  tower.  It  is 
tempting  to  attribute  this  period  of  low  z0  values  to  the  ocean,  particularly  given 
the  high  wind  speed,  but  if  this  were  true,  then  when  the  winds  were  from  the 
east  and  even  stronger  a  few  hours  earlier  (about  midnight),  z0  should  have  also 

been  low,  since  in  this  direction  the  tower  is  only  about  1 .5  km  from  the  shoreline. 
An  additional  factor  is  the  level  2  humidity  sensor  failure  during  this  time,  which 
results  in  an  unrealistic  stability  profile  and  has  a  minor  effect  on  the  magnitude 
(though  not  the  trend)  of  z0.  The  effect  of  backing  winds  due  to  increased 

frictional  drag  over  land  was  considered,  but  any  correction  for  this  effect  would 
require  that  the  120°-140°  flow  observed  at  the  tower  corresponding  to  the  low 
z0  values  in  the  early  morning  hours  of  day  232  would  have  a  greater  southerly — 
not  easterly — component,  hence  an  even  longer  fetch  over  land.  In  this  case,  a 


61 


SHF  (Wm 


more  reasonable  conclusion  is  that  the  low  z0  values  in  the  early  morning  hours 

of  day  232  are  probably  a  result  of  flow  over  the  cleared  land  around  the  Mars 
Observer  Launch  Facility. 
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Figure  38.  Same  as  Figure  6,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  temperature  levels  2 
and  3  of  Tower  01 10  on  18-21  August  2008. 
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Figure  39.  Same  as  Figure  7,  except  observed  conditions  and  calculated 
quantities  are  from  the  southeast  sensors  at  wind  levels  3  and  4  of 
Tower  01 10  on  18-21  August  2008. 

C.  ROUGHNESS  LENGTH  FROM  DIFFERENT  OBSERVATION  LEVELS 

Figures  40  and  41  are  time  series  of  wind  conditions  and  roughness 
lengths  calculated  from  observations  at  various  combinations  of  levels  of  Tower 
0002  on  23-25  September  2008.  Though  this  period  is  not  included  in  the  cases 
discussed  earlier  in  this  chapter,  the  variation  in  z0  shown  here  demonstrates  the 
multi-layer  structure  of  the  surface  layer  in  the  region. 

Figure  40  shows  the  wind  direction  and  wind  speed  from  the  northwest 
sensors  at  the  lowest  two  wind  levels  of  Tower  0002  on  23-25  September  2008. 
During  this  two-day  period,  the  wind  direction  was  from  the  northeast  with  speeds 
of  5-8  ms'1.  Figure  41  shows  the  roughness  length  calculated  from  different 
combinations  of  wind  level  and  temperature  level  pairs.  In  the  legend,  SE12/12 
indicates  that  z0  is  calculated  for  the  southeast  sensors  using  observations  at 
wind  levels  1  and  2  (the  digits  before  the  solidus)  and  temperature  levels  1  and  2. 
Three  groups  of  roughness  length  variations  are  observed:  z0  from  the  lowest 
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two  levels,  the  middle  levels,  and  the  upper  levels.  Values  of  z0  for  the  lowest 

layer  (blue)  are  on  average  the  smallest,  with  a  mean  roughness  length  in  this 
layer  of  about  0.5  m.  Roughness  length  from  the  top  layer,  between  44.2  and 
62.2  m  (magenta),  has  values  of  over  4  m  and  a  large  temporal  variation.  It  is 
likely  that  these  two  levels  are  usually  above  the  surface  layer,  and  thus  z0 

calculated  from  Monin-Obukhov  theory  is  not  valid.  Results  from  the  various 
combinations  of  the  middle  levels  are  consistent  with  each  other  but  differences 
are  noted  between  the  middle  levels  and  the  lowest  layer.  These  results  suggest 
that  the  lowest  two  levels  are  probably  within  the  IBL  induced  by  roughness 
elements  in  the  immediate  vicinity,  while  the  middle  levels  sense  the  boundary 
layer  advected  by  the  northeast  winds.  The  relatively  small  roughness  length 
calculated  for  the  lowest  layer  is  consistent  with  the  smooth  surface  along  the 
northeast  branch  of  the  "X."  At  higher  levels,  sensors  observed  turbulent  eddies 
induced  by  trees  and  shrubs  farther  to  the  northeast. 


Eastern  Standard  Time  (days) 

Figure  40.  Temporal  variation  of  observed  (a)  wind  speed  and  (b)  wind 
direction  from  the  northwest  sensors  at  wind  levels  1  and  2  of 
Tower  0002  on  23-25  September  2008. 
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Eastern  Standard  Time  (days) 

Figure  41.  Temporal  variation  of  calculated  roughness  lengths  at  various 
combinations  of  level  pairs  from  (a)  the  southeast  sensors  and  (b) 
the  northwest  sensors  of  Tower  0002  on  23-25  September  2008. 


During  most  of  the  daylight  hours  of  day  267,  when  wind  speeds  were  at 
their  highest,  z0  values  from  all  layer  combinations  are  relatively  uniform  in 

magnitude  and  suggest  that  flow  during  this  period  did  not  exhibit  multiple  IBLs. 
This  period  corresponds  to  a  wind  direction  roughly  aligned  with  the  northeast 
branch  of  the  “X”  at  this  tower.  During  periods  in  which  the  wind  profile  has  a 
multiple  IBL  structure,  Monin-Obukhov  theory  must  be  applied  with  caution  in 
order  to  avoid  crossing  from  one  layer  to  another.  Doing  so  would  invalidate  the 
“constant-flux”  approximation,  since  each  IBL  is  in  effect  a  separate  surface 
layer.  A  method  of  correction  for  this  issue  was  not  applied  to  the  calculations  in 
this  study,  but  is  worthy  of  consideration  in  future  work. 

Moreover,  Figure  41  shows  the  temporal  variation  in  z0  from  both  the 

northwest  and  southeast  sensors.  During  this  period  the  wind  direction  is  from 
the  northeast.  Therefore,  flow  distortion  caused  by  the  tower  does  not 
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significantly  influence  the  calculation  of  z0  from  sensors  on  either  side  of  the 

tower.  Hence,  the  calculated  results  from  both  sets  of  sensors  are  very 
consistent. 
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V.  CONCLUSION 


A.  SUMMARY  OF  RESULTS 

In  this  study,  near-surface  observations  are  examined  to  determine 
relationships  between  surface  characteristics  and  the  temporal  response  of 
roughness  length  z0  from  various  altitudes.  Although  the  WINDS  sensors  are 

capable  of  measuring  and  recording  observations  at  higher  accuracy,  relatively 
imprecise  measurements  were  available  for  this  study — particularly  the  wind 
speed  data — that  were  given  in  integer  knots.  Nevertheless,  after  averaging  the 
observations  over  20  minutes  to  better  represent  mean  conditions,  reasonable 
values  for  surface  fluxes  and  roughness  lengths  were  calculated  using  the 
WINDS  data  for  various  stability  conditions  when  wind  speed  exceeded  2  ms'1. 
In  low  wind  conditions,  the  coarse  precision  of  the  data  together  with  inherent 
uncertainties  in  Monin-Obukhov  similarity  theory  for  free  convective  conditions 
resulted  in  unphysically  large  values  and  temporal  variations  in  the  calculated 
roughness  lengths.  To  mitigate  these  effects,  observations  taken  at  low  wind 
speed,  which  occurred  most  frequently  in  the  nocturnal  stable  surface  layer,  were 
excluded  from  the  calculations. 

Two  periods  suitable  for  the  intended  analysis  were  identified,  coastal 
wind  circulations  on  1-4  June  2008  and  Tropical  Storm  Fay  on  18-21  August 
2008.  In  the  coastal  wind  case,  the  sudden  change  in  wind  direction  following 
the  passage  of  SBFs  provided  an  ideal  scenario  to  evaluate  the  variation  of 
roughness  length  with  wind  direction.  During  the  tropical  storm,  high  wind 
conditions  minimized  the  relative  error  due  to  the  precision  of  the  wind  speed 
measurement,  thus  increasing  confidence  in  the  calculated  surface-layer  fluxes 
and  roughness  lengths. 

Results  revealed  significant  variation  in  z0  consistent  with  upstream 

roughness  elements  to  a  distance  of  about  100  m  to  2  km,  which  varied  with  wind 
speed  and  observation  height.  Lower  observation  levels  of  the  tower  were  often 
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located  within  IBLs  induced  by  the  surface  in  the  immediate  vicinity  of  the  tower 
that  limited  the  upwind  distance  roughness  characteristics  could  be  sensed,  even 
in  stronger  winds. 

Ultimately,  the  goal  of  this  and  similar  research  is  to  develop  improved 
methods  to  incorporate  the  dependence  of  turbulent  fluxes  upon  the  ABL  wind 
field  in  mesoscale  numerical  weather  models.  To  illustrate  the  potential 
implications  of  this  dependence,  Figure  42  compares  observed  and  calculated 
data  from  the  WINDS  to  the  Coupled  Ocean/Atmosphere  Mesoscale  Prediction 
System  (COAMPS)  model  fields  for  selected  times  during  Tropical  Storm  Fay. 


Figure  42.  Contoured  roughness  length  from  (a)  COAMPS  and  16.5-m  wind 
observations  and  calculated  z0  from  the  WINDS  at  (b)  0000  EST 

and  (c)  1800  EST  on  20  August  2008.  COAMPS  data  courtesy  of 
the  Naval  Research  Laboratory;  image  ©2010  Google. 
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In  Figure  42,  satellite  imagery  of  the  Cape  Canaveral  region  is  overlain  by 
interpolated  z0  values  contoured  at  0.2  m  intervals.  As  previously  noted,  z0  is  a 
fixed  quantity  at  each  grid  point  in  mesoscale  models,  including  COAMPS. 
Figure  42a  is  the  constant  COAMPS  z0  field.  Across  this  region,  COAMPS  z0 

values  are  uniformly  low,  ranging  from  about  0.2-0. 6  m  (see  color  scale). 
Roughness  height  over  the  surface  of  the  Banana  River  is  assigned  a  suitably 
low  value,  within  the  limit  of  model  resolution.  At  the  resolution  of  COAMPS, 
subgrid-scale  surface  heterogeneities  such  as  those  considered  in  this  study  are 
not  discernible.  Figure  42b  is  a  contour  plot  of  z0  (using  observations  from  the 

lowest  two  tower  levels)  at  about  0000  EST  on  20  August  2008,  when  the  center 
of  Fay  was  located  southwest  of  CCAFS.  Arrows  indicate  the  speed  and 
direction  of  the  16.5-m  wind.  Figure  42c  is  a  contour  plot  of  winds  and  z0  from 
18  hours  later,  when  Fay  had  moved  northeast  of  CCAFS  and  winds  were 
offshore.  The  earlier  image  shows  that  for  WINDS  locations  near  the  coast,  z0 

values  were  comparable  to  the  fixed  COAMPS  z0  field  in  the  strong  onshore  flow 
conditions.  In  the  later  image,  with  a  strong  offshore  flow  component  and  about 
the  same  wind  speeds,  areas  of  higher  z0  values  are  present.  An  area  in  the 

center  of  both  plots  was  excluded  in  order  to  avoid  interpolating  across  a  large 
area  devoid  of  observations. 

Compared  to  the  COAMPS  gridded  data  field,  relatively  few  observation 
sites  were  used  to  construct  these  contour  plots;  however,  values  near  the 
higher-density  tower  locations  carry  a  high  degree  of  confidence,  such  as  those 
along  the  immediate  coastline  and  in  the  cluster  of  towers  located  at  the  top  of 
the  images.  The  observed  temporal  and  spatial  variation  of  roughness  length  is 
not  represented  by  COAMPS,  which  depends  upon  land  surface  type  only  and  is 
independent  of  wind  flow.  Relating  the  roughness  length  to  the  surface-layer 
dynamics  would  involve  a  complex  blending  of  the  effects  of  various  IBLs  over 
heterogeneous  regions,  and  this  relationship  would  be  necessarily  dependent 
upon  wind  speed  and  direction. 


69 


B.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH  USING  THE  WINDS 

OBSERVATIONS 

Wind  speed  observations  with  a  precision  of  about  0.05  ms'1  is  measured 
and  recorded  by  the  WINDS,  but  was  not  available  for  this  study.  Use  of  this 
high-precision  data  would  have  greatly  decreased  the  number  of  exclusions  due 
to  insufficient  speed  shear  between  wind  layers  and  permitted  more  numerous 
and  accurate  calculations  of  roughness  length  and  surface-layer  fluxes  for  all 
cases.  For  future  work,  it  is  recommended  that  this  high-precision  data  be 
acquired.  At  this  level  of  precision,  more  careful  consideration  of  flow  distortion 
due  to  the  tower  structure  would  be  required.  The  potential  for  this  distortion  to 
affect  the  calculations  was  recognized  in  this  study,  but  no  algorithm  for  dealing 
with  its  effects  was  developed.  Research  into  the  effects  of  z0  variability  on  the 

exchange  and  drag  coefficients  and  how  these  effects  should  be  simulated 
efficiently  in  high-resolution  numerical  models  would  also  be  a  beneficial  study, 
as  would  the  effect  of  clouds  and  precipitation  on  z0 . 
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APPENDIX.  WINDS  TOWER  LOCATIONS  AND  INSTRUMENT 

COMPLEMENTS 


Table  1.  Locations  and  instrumentation  heights  for  the  four  launch  critical  WINDS 
towers  at  CCAFS  used  in  this  study.  After  CSR  (2006). 


TOWER  # 

LOCATION 

INSTRUMI 

ENTATION  HEIGH 

T(m)1 

Wind 

Temperature 

RH 

62.2 

62.2 

62.2 

44.2 

- 

- 

0002 

28°  26'  39"  N 

27.4 

- 

- 

80°  33'  44"  W 

16.5 

16.5 

16.5 

3.7 

- 

- 

- 

1.8 

1.8 

62.2 

62.2 

62.2 

0006 

28°  30'  47"  N 

49.4 

- 

- 

80°  33'  41"  W 

16.5 

16.5 

16.5 

3.7 

- 

- 

- 

1.8 

1.8 

62.2 

62.2 

62.2 

0110 

28°  34'  1 1  N 

80°  35'  12"  W 

49.4 

16.5 

3.7 

16.5 

16.5 

- 

1.8 

1.8 

150.0 

150.0 

150.0 

120.1 

- 

- 

89.9 

- 

- 

03 132 

28°  37'  32"  N 

62.2 

62.2 

62.2 

80°  39'  26"  W 

49.4 

- 

- 

16.5 

16.5 

16.5 

3.7 

- 

- 

- 

1.8 

1.8 

1  Launch  critical  towers  support  dual  instrumentation  packages  at  each  level, 
aligned  northwest  to  southeast  on  Towers  0002,  0006,  and  0110;  and  northeast 
to  southwest  on  Tower  0313. 


2  Tower  0313  also  houses  redundant  Vaisala  PTB220  Series  barometric 
pressure  sensors  at  1 .8  m  AGL  (4.3  m  MSL). 


71 


Table  2.  Locations  and  instrumentation  heights  for  the  14  safety  critical  WINDS 
towers  at  CCAFS  used  in  this  study.  After  CSR  (2006). 
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